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- -S3 

where /i is the interval of differencing. Let fax) be a polynomial of the nth 
degree in x taking the same values, as y corresponding to x =*. j.e;, 
^/-represents the continuous function y - /(xj such-that /fx^ “ r = 0„ 1, 2, 

m and at alkpther points f(x) = fax) + R(x) where R(x) is called the error term 
(Remainder term), of-the interpolation formula. Ignoring the error term iet us 
assume / . 

f(x) ~ fax) ~ a Q + (x x^ + (x - x^j (x — Xjj + ... 

+ (x “ x Q j (x - x } )..: (x — x 

the constants a Q , a y a v ... a n ca'n ; be determine as follows. 

Putting x = x Q in (3) we get 

f< x o> - r a 0 .. ' 

putting x = Xj in (3) we get 

/<V ~= a o + a i (x i ~ V - Vo + a J A 

y, = y 0 + a i h 

1 . A ■“ ' 

Putting x = x^ in (3) we get 

f(x 2 ) = (/x'jr-.J = n 0 + - r/CXj - a',} 



YX?! 



•-y 


/ 0 + -fa 2h ) +a i( 2h i h ) 

h 

Vi = y 0 +2 (vi - yo) + ‘ I 2( 2; ' 2 ) 


y 2 - z y\ + y B 


y« 


- 2 2ft 2 2 !ft 2 

Similarly by putting x = x y x = x 4 , ..., r ~ x^ in (3) we get 

yo 


X". 

Ci> 

Her 


P 


^ 3 y« ^y n 


3 3 ' 4 4 

3!fc 4!fc 


y 0 ^ yo 

.... a — -- 

nlh n 


putting the values of a Q , a y ... a* in (3) we get 

f(x) « 0(x) = y 0 + — (x ~x 0 ) + — y(i - jr 0 )(x - .r , ) 
U 2\h 

+ 4^^ “ X ° )( X ~ X ')( r “ X 2 )+- 

3 !n 
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m 


.f 


- = 1.6 


We have = 12500, h = -iQ.arid .x =* 12516 

/jr-jc 6 . 12516-1-2510 
v . . h . . JO. 

from Nevvton's forward mterppl at ion formula-/ - , 

/( x )=yo. : t^yo + ~~~ ; 1 ^ zl2 yo + - 

12516» fl fc8f&4:*- -1:6 -X>0:0447 4- 0 
. " ; = -111.8034 ’+' 0.07152 

= 111.87492 >■ 


1.6 =111:87492. 


Example 5 : Evaluate y = e 2 * tor x = 0.05. using the following ta&fe 


y 0.00 , 

0.10 

0.20 

0.30 

0.40 

y = t.OOO 

’ 1.2214 

4/4918 

1.82-21 

2.255/ 

Solution : The difference table is 


— 


N- 

11 

Ay 

A 2 y 

4 3 y 

dV 


x “ x o _ 0-05 -0.00 


0-1 


—05 


Using Newton's forward formula 

/ \ i^(u — l) , u(u~— ljfu-2) .. 

/(*)= y» ! "4y 0 *— - A y 0 +~—“r 


u{u — l¥u —2 Vu —3) 

+ ^^ d 4 y +... 

4! 0 


3V 


ss 


3£ 

B 

1 

2£c 
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-?• 


0.000 

1.0000 
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Cr 0 i 


0.2214 



.1 


i 

0.10 

1.2214 


' 0.0490 



" 




0.2704 


.,0.0109 



s 

0 20 

1 491« 


0.0599 

0 0023 



1 



0.3303 


0 0132 


■* 

1 

0.30 

1-8221 


0-0731 







0.4034 



, “ 



0.40 

2.2255 



— 1 












0.0.0, x 

= 0.05, 

ft = 0.1. 
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1S009+O.5X 02214-fr f—^ v —- (0.«490)? • . 

' ,i ;:-2... ‘ 

= 'B.OSS “ 0.006125 + 0.000681 0:000090 

==*I«!6£... . 

1.052. ‘ 

S<iamgfe-6c -.iTSfe*:. raises of sin* are gixen beteMv.for drffe^&^:vShi'es;.€ifi ' : §mdi 
b&e r v;adke‘;f&.sam 33° :- . 


•£ 


30° 

. 55° : 40°; 454 

5CT 

14 

y = stiur 

0.5000 

0.5736 0.6428 0.7071 

0.7660 

ri 

Safutzoib .^ot.. — 32** ts very 

near to the starting value x ' ; 



rbr^vard- irt^R^afergiiDTs formula the difference table is 



X 

q =smx 

Ay A 1 y A*y 

dV 


30° 

0.5000 

0.0736 



0 35* 1 

0.5736 

-0.0044 





0.0692 -0.005 


=' 


0.6428 

-0.0049 

0 




0.0643 ^-0.005 



-£5° 

0.707] 

-0.0054 


* 



0.0589 


s 

50° 

0.7660 

■ 



_ x ~ *o _ 32 ° -^0° 


= 0.4 


A .5 • 

SXfcfcaveg^ — - 0.0736, d^ = -0.G044, d?y 0 ■=> ^0.05^- : '--^. : - 

psa&jag. t&esses isi We wEen's forwardd.aterpoiatidn^rd^roXr^'v^giE^ 

~ T ) . “( K - l )( U ~ 2 ) f :. ■ . 

-2J y 0 +-;-— Ay ! * ■ 


■*%=»&+* 


2! 


3! 


»Xp2.9|j=€iSSJa& + 0.4 X 0.0736 4- i — iX^ 4 • *) (^0.0044),- -*:.«•.* / 

— aseoo + 0.02944 + 0.900528- e.ooe.32 = 0.52993& - #29^. 

0.^*8 


i 


fe Sfcr 


It 

£ 

I 
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JExample^^/fn an examination the number q(. candidates who obtained marks 

between certain.limits wer e as follows : 

Marks 30-40 40-50 50^60 \ '6&—7Q .70^0 J: - . 

No. of Students 31 ' : 42 5*1' ' *3 5 ' ; 31 ' 

Find the number®?/ candidates whose scores lie hetween 45 and 50. 

Solution c. First, of .all we constructs xrumuiative'frequency table for tire given 
data. . . . . ■- ;■ - • ; 

Upper limits of the cjass intervals 40 -50.. .. 60... ,70 ."..-3D 
Cumulative frequency 31 73. 124 159 1 90 

The difference table is < 


X 

y 

Ay A 7 y 

4‘\y d‘y 

marks 

cumulative 



frequencies 


40 

-31 

42 


50 

73 

9 ■ ’ 




51 

-25 

60 

l 24 

-16 

37 



35 

1 2 

70 

159 

-4 




31 


80 

^190 



>, h = 

10 







u 

_ - r . - *« = 

* h 

45 - 40 - 

10 

0.5 


andy 0 ^ 4y 0 = « = 9 - = -25, tf y# . = 37. 

From Ne wto n's forward interpolation formula 

’ 

./ \ . — l) u(u - l)(w _ 2) _ 

/O') =- v 0 + "4v 0 +■-— +-L—Ji- -±3 y fl 

• * 

u (» - Ob - 2 )(“ ~ 3 ) 


_ +- 


4! 




- /- /(45) =31':#(0.5)(42) + ( 0j5 X °- 5 ) x 9 


(0.5X0-5 - l)(6-5 - 2) / 
6 ' 


(-25) 

+ (u.5X0:5-l)(o.5-2)(0.5-3) ^ 


24 


: : /* 

7 fc 
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' — : S'jt2Sf-— IJ5625.-3A452:- 

=47J8373 s . 


-•The number. of. fee 

Rsm&er of ^^CTB-.wKo. s^e<J E*arfcs 

SSfefion : We c^fe^ differezrre iabife 


v-i-r- 

\> lit 

- - " _x 

' y 

*dy 

4S 


r “i 






M 

2 

~r 

0 

t 

2‘ 





2 


0 

ff 

3 

x 

4 

. 2 



4 

5 





We have = f,ft = l„y 0 - -1, = 0, 4 2 y 0 = 2, 

x - r 


* 




(?Eom New-f&n's forward interpolation w e get . " 

y ~ f{ x ) ~ yo + “Ay 0 + -^—^-A 2 y 0 +. 

2| u 

^ ^ r ) = • + (x — l) o + —— Lull .2 

2 • 

ftx) — X 2 - 3x + 1. ~ 

SUV ?T “? *? POiy nC>roiak ° f hi S»- degree ^ik dat^fc*. 

taBt8,B WS ' 8fte P° l y«°«»ial of least degree^fc^^te d^f.v 

•.«§£»»« ilip^fei:W3Fd mrerf3Q!a(jBB fsrasfa^a . ; T. 

^^.Sb^- i^s^a^n formula caanot-be^ed^^fe^^ a 

f * taMe of Yalttes - ** ^p*h^..^w 3i&i L 
ST?fT 1 VV ^ £ ° n ' <Gree ° ry 

Be derived 25 , fbM&ws: " ■■ . 

» S .-JMta-faMtaa „HUk „te .he values v coa.esp.iute,. 

he equaJiy spaced! wife h as the interval of differencing,..j, Q;r : - 
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yi 

■n 


o 


ft 




**-3 = 1 " C-.* ?» 


(uh) + l I -€■ I# ; ; 

X — X^ 7 =. <» (* — Xji = ( !i.+ ^ — 1> it ^ : - 

AV^e-'etjxiaifozi- 66> may be wriffen. as-\ ' ;" . -..t^; ' .- : ::; •■.. 

&% pg - * T li - . + t£ -2) ... 

■- = + -^ + ^-~V 2 y ^ ' ' ., c. 

^v_.__ r ■ / . p 2 f • • *----■ 

c + c t ,> ~ 1 ) » % • 


The above foymula is. krtovw as^ Newton's backward ircierpola• 

/ .... /,.. * 

BxaHrjde 9 : foUowing data gives the meeting point of a>xt allay;.o5'Jea%ap^f; 

zinc, where, f is the temperature- in degrees' e and P is the perceafe'ge af : feadr; in 
the alloy. , 


& 

40 

5CS ■ 

60 

70 

80 

90 

& 

ist§h 

aiu 

226 

250 

276. 

- 304 


Find the melting;- point-of the alloy Gontaining 84 percent lead'. 

Solutioirr: The value of 84 rs near \t$ie end of the table therefore we. ose ^h.e 
Mev/tons hae^v/ard interzrolatiou foximda. The difference table is 


F 

* 


V s V 

F* T 5 

40 


2.6 



50 

2.(54 


2 




22 

,0 


60 

_226 


2 

-0 



2*? 

0 

0 

70 

250 


2 

0 



26 

0 


SO 

276 

28, 

2 


90 

35.4 





We have x* x = 84v-&=l$>.1* — y rt =304, V{* = Yy^ = 28, ^7 2 y ^ — 2', anc£-'' . 

- . ^ = ^=^.=0, 


i-x . 84-9Q 




18 


= -4).6-: 


From Newton hac^waTtl fomurla- 


id-v -i~ X) _» 

f £84} = i + uVi -f—- V i 

‘ - - 2 


3 r 


1 


I 

1 


m 




1 

i 

i 


w. 
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Numciicdf ArrdljtiS: ■ 


lied -the fucEfetkm. ftx > be coatnwoos astd: possess continuous cterivatiYes.of all 
otrdds wilii m the interval [i fl , "rrjf an^ Jef.^fcxJ denote Ahe .interpolating 
^srtynome fehe auxiliary iunctton'FH) as given below^ • 

• U i ; . - ' . .I 

■* - •; - i* ^oA* ^j) -u x «) 



- Wit) satisfies all the conditions of Rogers Theorem in each of/theiV^^aitetyafe-- 
fix ;x x ), ( x y A x n -i' .* JK H^nceT^fi vanish^s at least once in each of fcfie.-»jr&- 

mtervafs. Therefore f vanishes at least (n + 1) times in (x 0/ x^V.^' r (£}vaidshes 
at least n finies m the interval f* x ), F a ^(t) vanishes at least once in? 

■■ «■ « ' . * -r- /* 

gx e , x,^ say at £, where * 0 < £ < x n - 

The expression 0 ft — x y > (f — xj is a^olynocniai of degree .fw + 1 ) in i 

and the coefficient of t — 1- * 

_\ The (n -f ljth derivative of the polynomial is fn 4- 1>! ~ 

I- + 0 ! 


•■• f' +l (^) = /‘ r+, (^-{/( 1 )-^)}; 


(x -x 0 X* -*.) 

=> /(*) -<*h ) - {" J|jh - -'c/i* - r 1) ■ - (* - r J 

111 Rfx; denotes the error irv the formula then R(xJ = /(*)." $(*) 

K h)=(* -jr *X* - * >) - ( x - *,) 


— O' 


But x — x, Q = tch 


— x {~ f« 


(n+l): 

u ~ UJ 

of differencing therefore we can write 
- T- EiTor 


Tlfc, ... fx - xj — (u — nj/i where./r is the lfitfervai 
4 




Using the relation D ~ — A 
h 

we get 


h H * x 








A 


) * 
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{$ 

S 


; r ~t 

8 


gf 


1 




n + X* ' 


The error in theTorward interpolation formula is - - - :■ ‘ 


-9( u - 2 K( u_n ) 

(»+ O’ 


Similarly by taking the auxiliary function FU) in the form 

and proceeding as above we .get the errot in the Neiwton 'baefcwarcl 
interpolation formula as ' - \ 


r(x) = +!)...(« + n) where uh =x.~x^Tj 


Example 11 : Use Newton's forward interpolation formula and find.the value 
of sin 52° from the following data. Estimate the error. 


I 

-15° 

50° 

55° 60° 

y = 

sinx 0.70 71 

0.7660 

0.8192 0.8660 

Solution : The difference table is 

r 

V = sin.r 

4y 

<4 7 y d 3 y 

45° 

- ' 0-7071- 

0.0589 


50 ,J 

1J.766U 


-0.0057 



0.0532 

—0.0007 

55° 

0.8192 

0.0468 

-0.0064. 

60° 

0.8660 



We have x Q = 45°, x, = 

; 52°, y 0 = 

0-7071 

, Ay 0 ~ 0.0589, £ 

4 3 y 0 = =0.0007 , 





■ X - x 0 

52° - 

45 ° -14 


u = 

A 

5° 

-i- 1 - 


From Newton's formula^ \ . 

ulu l) ^ — l)(u -2) 

I/=U +I«A' +-zfy + — ^y 0 + - 


2 \ 


3\ 


* : 

■'i 


1 

£* 
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& 


3 

•v ? 

^ I 




*4 I 


I 


rx 0 : 0589 ^ *). / 

■ x f Q.0Q57): 

+ ih^±^rA ~^ f .. 

6 ~~ 




*"*■ sin 52° = &.78Si0Q3-2- 
Error = 

etldcig: » ; =■ 2 we gel 




E> ror =jfe r 

' /X 


3f , 




- -T)f?,4 - 2) / ^ 

■6 


= 0-0000392 . 


11 


Flnd missi ^^ t^foir owing tafe!e 


fix) . 2 


3 4. 


* ssr “* ^«««« o f „. 

year <7T 


* fi ° d ""’ 

E 5 E 1 ! ggta 

*7‘” d ““"**■ »■ : 

* *• « ; ^-■ 


—15 20 .2-,s 


templet* i^e- taoble 


* ' ^nrf/n.p 


t8 


32 





'i: 


M -i 


m 


I r 


M 
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~~ .z/j — . 5^ i tg == 52Q ,Uy = 3S5 show that 465. -> 

m fa an exii ntina fipn. fcjti^ .nuuiber:. of candidates who Secured otarks between 
certain limit .were.as idHrmS :. . •'' ' 4 " een 




. N^fol candidates- 41 


ft 19' W ~ 39 40 ~ J59 60 —79 SO — $9 


63 


■65 


50 


• 17 


' cindffdaW W fe^e inar^>Ve -less-th^V 

29- Given the fo.Hcwtijg score diatribuffaTi of statistics 
- Mark s ■■-■%. 3ft --■«y--«>. ^5b-5)--60;; 6 ^I.W 


No'oi Stud:en.ts -52 


21 y 74 


Find 


(i) the number of s students" wMo secured below 35 
(it) the number,pf students who secured above 65 
<iy) the number of students who secured between 35^45 


%ssrs 


»- 17 2. 6-60 3-?8 4- 23.5, 1.4.25 5. fO) = .0.123 f(6) = 0 090 fi -t « 

2 aS9 9 8.5281,6^04 9.10.963 1 0. 2,0959 11. 0.46889 12 27992 

®./ + r + 41 14.4*3+, 15! 2*3 + 3a -1 16. * 2 + * + , 17 ,3 _ £ ‘* 

39. 197 20. (i) 26 (ii) 7 (Hi) 46 ; 2 + “ 


^o-bU 4ijL. JuSfdnCg,. in wftUKCa^ vv\Cf<LS a| 

M-o >.?-+? U-» t_ ItL — dt „ „o ., . 


Vl‘i ' y * h OY'+i-'Ori ' f ‘ 0 T Cf?\,K+s 

° TJ (J 

abovd. Gfxxftii t'iueflo.tj. 

. 2.&0 


‘ . '/• ' 1 :—~ --—v*. 

JsJt^U- loo iSb. . %o.o • sufQjj*'S.00 c f oo |. _ 

( c ’-6^ .12-03 fS'-O^. 

(_ A\ - . . - ix*' - t- 

' —--- —__ _ _ . £*■/ 


I'Cvul “tbe vet 

v 

(rj -x r 2-( $ 

w j-mo .ft- 


a 


- 


3f. 

il 


SI 






I .i 


4 


g 

& 


I 

£3 


& 

i 

B3 

I 
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-# 


we get __ 

no-**#*-,,* 

.(5-6)(5-9X5-ll) _ > 

^(to-5Xie-6Xio^. y ^fl?-**;**-*. . 

- ,p-5)(9-6){9-n) < > 


_ 13 35 16 v 

= 2 :-(- -- 1 -- 

3 3 3 3 






$ 

l 

•i 


B . Find by Lagrange's formula the TaSernrdafi^rr ■ 

corresponds to the following data P P ° ynomial urhhrB 

* £ * -1 

fix) 


* Fin<J thC P Q] y nom *« d egree three relevant to the following data 


* -i ~o ~ 

J . 2 

. /<*•; l i 

l -3 .. 

Find the polynomial of the least degree 
a l the given points 

which.attains 

* -V -2 1 

2 4 

fix) 25 -8 

15 —25 

Compute [(0.4) using the table 


* 0-3 0.5 

0.6 

fix) 0.61 0.69 

0.72 


<» of a c 

thej>aJes for the years (a) 1986 (b) T992 


Vear J985 
, “ Sa Ies 40 


‘ _ (in thousands) 

Compute sin39? form the table 


m7 ms 
48 


43 


T99T 19 93 
57 


52 


ia 


si r>jc° O 


20 - 


3(f 


40 


— ^ F ^ nd Iog^5:I5 from the table 


l .1-736 0.3423? 0:5000 ~ 0_642fl 


St 


91 


x,. 5.1 

5.2 

5-3. 5.-4 ^ - 5.5 

,0 S, 0 * 0-7076 

0.7P60 

0.7243 0.7324 0.7404 

s interpolation formula arid find 1(0 35 ) 

x 00 

0.1 

G-2 0-3 o 4 

f(xX. 1.0000 

1.3 052 

1-2214- 1.3491 i aq\u 

—-----y i.o 


Use Lagrange's formula and compute 


v: 


£ 

* 


s 


S 

-6 
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5 'lr 


a 

m 


£2 


y 


as'; 


Inici-fiolation iurtk unequal -Intervals 

9. Use- tograng;e-s formula arid- compute 


X . 

0.20 

0.22. 

0.2.4 0.26 1 

0.25 . 0.30 

fU) 

- 1.659.6' 

1.6698 

V.6804 1.69:12 

1.70Z4.-1-^139. 


» . » . - . ; ; 

10. Use Lagrange's interpolation formula to :find the value b£ fix) for jc — 0 




r~".-v : 

' -2- 

- 2;\ 

4; 

fix)] . 

. -9 - 

- 11 

•69 


ARSwers ;V : . 

1 .x*-3i +5 2,,-or? + Jf-+.V -3. X 2 --.10x +1 ■ 4y;^.65 L 5. (a)4IJE2 ^^54.46 : 

6. 0.6293 7. 0:7118 8. 1.4191 9.1.6751 10. 1 ■ _ 


S 2 < 

*Ss 


© © 
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Bsearmpfe : The follow tug table gives the value of^ the: elliptical ; integral 

certain of $ . Fired. ;the values;.of $ if T(<£) : ~ * 8:3&&7 


- 

-’ll 


21 ° 


23° 


25° 


Ff/> 03706 6.4068 0.4433 


Sfrlulicw* : We. have 

# = 21?;*;^ ^.23°, <p. z < 25V £ 4= &3S87, f fl - O.3706;VF, ^mUSS 'and /* 2 =p;4433L - 

0/5ing the inverse interpolation formula we cari -write 

~(f-^)(f-4) V 

- K-^)k ° ft-jife-*)* &-*)&-4)**' 


‘"■s 

.£71 :- 
T 

l! 

v- -i 


(0.3887 - 0.406g)(0.3887 - 0.4433) 


x 21 


(0.3706 - 0.4068)(03706 - 0.4433) 

( (0.3887 - 0.3706)(0-3S87—04433) ^ 

(0.4068 - - 0.3706)(0.4068 - 0.4433) 


x 25 


(0.3887 - 0.3706X03887 - 04068) 
• (0.4433 - 0.3706)(0.4433 - 0/ 

= 7.864 + 17,20 - 3.087 
= 21:999 22 
.'. $ = 22°. 


Saastplg , : Find the value of x when y = 0.3 by."-applying Lagrange's forn>ufa 
inversely .V.,7. •- 


X 

0.4 

0.6 

S.8 

y 

03683 

0.3532 

0.2S97 




Salution : Frefti-Lagrange's inverse interpolafeson-formula we get 

(y-yi)(y-y z ) "7 (y-yoK^-y?) (y-y«X y~yi ) 

(yo-yi)(yo “Fz) ° (y>-yoXyi-yJ * (yz-yoX^"^) 


' SotsstjJufing-x p = 0-4, Xj — 0-6, Jr 2 = G.&, y G = 03683^^^ == 0:3332, y 2 — 0..2&99 in the 
r£> -.v&ove formula, we get ~ 


fi 

I 


i 


5-1 


u 

U 


1 
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?.ct 

M 

XI 


fer 

rSSj...- 

.%3fei ■ 
JgSV. 

. sr -' 


3 


3 

|| 

v5 


3£ 


Interpol a Uoffivith unequal Intervals ' 

■ (03-o:to2)(o;3- ' 


V* Uiv 


(0.3683- 03332)(0.3683 - 02897) 

- (0.3 —0.3683)(0.3—02897) 

(0.3332 - 0^3683)(0.3332 -02897) 
(063 - 0. 3683J(0.3 — 0.3332) 
(02897- 0.3683)(02897 -0.3332) 
= 0257358-. 


X (0.4) 

X (0.6) 
X (0,8) 


(J^ample : The following table gives the'values of the probability'.intggrb! 

2 r _ x 2 - V-i>5. * 

y = —-== € . rfx corresponding to certain values of^c. For what valuer^iJf is. this 
JjV J . ■■ v . • 


integral equal to ^ ? 


X 

0.46 . 0.47* 

0.48 0.49 

1 r-*\ 

y = —], dx 

v/r o 

0.484655 0.4937452 

0-5027498 0.5116683 


Solution : Here x Q = 0.46, x } = 0.47, x 2 = 0.48, =^A8, _y Q = 0.484655,. 

y, = 0.4937452. y 2 = 0.5027498, y 3 = 0.51-1 
From Lagrange's inverse interpolation 

v ' if, 


(y-yiKy-yzXy-yj) , (y 

— - -:*n + 


-yj)i 


(yo-yiXyo'yzXyo-'ya) ° fa ^g^pS^^ifa -y*) ' 
(y-yoXy-yiXv-yj) * (y-4sXy~yiXy - y2) 


-x, + 


(yz yoXyi'yiXy-z-ya). .-fa-yoXys-yi-Xyaryz)-*. 

(0.5 0.493?452)(t?5 -0.5027498)(0.5-0-5116683). , 


rfe-X 0.46 


(0.4846555"— 0.4937452)(0.484655 - 0.502.7498)(0 48465SS - 05 H 6683) • 

" (05 - p.4846555)(0.5 - 0.5027498)(0.5 - 05g 16683), -; , • 

— (D.4937452-0.4846555)(0.49374S2- 05027498)(0:4932452-0^tli|583) 

(05 - 0.4846 555) (05 — 0.4937452)(05 - 05.116683) . 

(05027498 — 0.484655)(05027498 — 0.4937452)(05027498 —05116683) 

t (05 - 0-484655)(05 -0.4937452)(05- 05027498) ^ ^ 

(05116683 — 0.484655)(05116683 — 0.4937452)(05 — Q5027498) . 

= -0.0207787 -tO.l 57737 + 0.369928 - 0.0299495 = 0.476937.. 


X 0.47 


X 0.48 


SS 

iCC 
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Solution :••• 


i 

<jXx).r 


v* 

5$ ,-f 

zx ■ 

H*) 


Vji 

i 


d 

8 


o 

SJ 

i5-2 

-•• a*.). = 

tw 

&§ 

" /(f 


d 


& 


Example 1 

o 

d 

(i) y, 

gcf 



A 

■T 

00 y. 


M Smm# : 

where (x- x^x x„) 

an^ = 

ax 


= (x — x t Xx — xj(..Jx ~ X-J ; 

+ (x — XqXx. — Xjf.-Cx — X o Jf 

+ ... + (V - XpKx — x Y )„-(x — x^Xr - x r J )...(x — x) 
+ ... ‘+0r — x^x 






Solution : (i) y_ 5 / y_ 3 / y 3 and y 5 are gtven> therefore the values of'the 
are -5_—3, 3, and 5, y-j is to be obtained: By lagrange's. formula 


y x 


[x-(-3)](r-5)(x-5) ^ r [x-(-5)](r-rf.3)(j.-;5).: ^ 


5—{—3)^-5—3)(—5 5} “ |-3-{-3)j{-3va)(--3,-5) 

[3 - (,-5}|i-.(~3)|3-^ 


ry-3 


Taking i = 1, we,:ge^ 

■STvS ' 

(1+3){1-3Mt-5) 


?! = 


. (l^5^l.-3Xl-5). 

y-s + , ^ „w, „ ^ y-a 


(-5+ 3X~5~3X r ^ 5 ) {-3-f-^|-3-3)( -3 -5) 

, (i + 5)(H-3)(l-5) t 
(3 + 5)(3+3)(3- 5) (5-f 5)$5+ 3}{5-3) ?S 


•«g 

• ^rZ 


i 


il;2 . 

ample ~ Show that Lagran-gef.s. inlerpoiatipn ;f:oruwiLa ; caiv be writers fn 

5ke form - ' : f -V . ; . 

U :.h \ • 


:<\*c 

^1' 

iSf 


_« V- 


I £ 

I 1 


& 
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§ 
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rrt 

*3 




for maarimifu, or Fnimcruiro we j,av e 

■*' A jg^fr+cj. y ( £_ 

=* H$t c M4+ (c-p)f (*)■+ (a -r 6)/(,r)} 

._•_ >* - 

■ + (c- tt)f(b-) -h ( a . .. 

r / . 

Example J . : Given lo^ 654 = o 6' 4 « f *. 

tog w 66l' = 7.8202’ firvd Iog I0 656. ‘ ' ^5^659 = 2.8589, 

Solution : Here * 0 = 654, x =: 658 . y _ /> - 

1 " 2 -.. *3 ~ 661 : Aad fix^ “ 1 ^^ 656 . 

By Lagrange's formula we have 


/(*) = 7 


1 

>%: 

K** 

1 

[* -*.0 



1 

F« - v ,j(*o -*,} 

{- V 0 - - r 3 


^-^,-^](x 1 -x 3 

r 


- **•.« = X ^ I% . 

{654 - 658)(654 - 659X654 - 661) ^ l56 > 

—*■ . ^656 — 654)(656 - 659^656- 661) 

(658 - 654X658^^^1^ X (28182) 

, ( 656 '~ 6S9)(656 - 658)(65 6 - 661) ?- 

^ (659 - 654X659 - 658)(659^66l) X ^ ZSl89 ) 

, ( 656 ~ 6 ^4 )(656 - 658X65 6-VWrt 
•’"- f 661 ~ 654X661 - 658X661^659) X1 82 °^ 

= ^( ? .8156) + |{2.8,82)-i(ii 18 ^(2.g2®2> ^ 

= 0.6033 + 7.045 - 5.6378 + 0.8 057 
= 2.8170 . 


u- -||g 


->3 

.Si 


|j 

il 

li 




i 




& 




BE 


cF:' 
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Interpolation with-uncqital IniertraIs ; “ 

Example. Write down the Lagrange's-polynomial passing through.(* 0 , /„), 

3x* + x 4-1 - 

far,,/,) and (x V / 2 ). Hence express (* _ j)( x _ 2 ){, _ 3 ) as sum of-partial fractions. 

Solution : The Lagrange's polynomial through the points (x ? , /„), (*,,./,) and 
(v 2 , / 2 ) is givenhf . 

/ / (x-x,)(x-x 2 ) (x-x 0 )(x -x 2 ) x + (*-*o)(* . 

/ X (* 0 -*tX x o'-* 2 ) " / (*z-*o)(*2-*i) 

Consider the numerator 3x 2 ^-+l. --J? . 

Let 

/Yx.) = 3x 2 +x+l, ’ 

tabulating the values of /<*) for x - 1 , 2, 3 vve get 


X ■; 

1 

2 

-7 

.* 3 

f(x) ' 

5- 

15 

31 


Using Lagrange's formula, we get 

/(r) _ kzlXi^l X 5x is + LzlXiz3l x3 i 
’ (l — 2)(l —3) (2-l)(2-3) (3-l)(3-2) 

= -.(a - 2)(x - 3) - 15 (.t - l)(r -3) + 21 ( x - l)(x-2) 

2 2 _ . 

. + l - = £( i -2)(x-3)-15(x-1)(x-3) + ^-(x-1)(x-2). 

(x-l)(x-2)(x-3) 2 - 2 


Esercise •< 

1 . Show thSt the sum of coefficients of y. s' in the Lagrange's interpolation 
formula is unity, 

2, -Given u r y |)7 u } , and u 2 _ UsingJLagrange's formula show that 


/ fy 2 — l) x(x 2 — l) 

— -i-- / ^ u ct V/here x -r y = 1 - 


u — i/u n 4-xu, -4— 

x J o \ 3* 

3. If y 0 , y y 2 , y 6 are the consecutive terms of a series then prove that 
y 3 = 0-05(y 0 + y 6 ) — 0.3 (y 1 + y s ' + 0.75(y 2 + y 4 ). 


, «; ■ • 
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o 


u 


3 


B 




*■ tf . «eept-fc. Of the-seqtWe^ ^ v y .,,.y 9 be given, show.^S' 

varne or is -- ■'■.*■ 

* 56 t. y «- + y'd ~ 28 (- v 3 + y 7 ) + s (y, +yj-(y, + y 9 .) 

70 - .. . 


3. Th e, folio wing table is given . 


.X 

. 0 

1'2 S . 

fM 

2 

■ 3 .12 147 


show that the form of f(x) is i 3 + x 2 - x + 2. 

6. Using -Lagrange's interpolation formula, Express the functior? 
* 2 +x ~3v 


X 3 — 2x ~ —x + 2 


as. sum of partial fractions. 


x 2 +6x +1 


7. Express the function f V7 '—77--- -* c CTtl ^ ^r. ,. , 

(-v — l)[x + l)(j — 4][x ^6} aS SU1TI °* partial functions. 

S. Using Lagrange's formula show that _ _ 


On ■: 

x- 1 — 1 Ox +13 

0 i 

( v -2)(.v-3) ( 

GO / 


Oh ■ 

x 2 + 6 t + 1 

H .i 

( A ' r -4 K a ' +6 ) 


2 _^ 3 


■ -2 


41 


9. The-foftpmag values of the funetionj^) for values of * are^iven: 7(1) = 4. 
/(2) - 5,(0) = 5, f(8) = 4. Find the .Values of f(6> and abo the v®hse*©f x fiJ 
which /U) is maximum or miniimim. 




; W|.' 

1 I 

! £1 . 

: 'a 

6 0^-l) (x- 2) : 2{; 


• 2 4 

?. “ ‘ + 


41 


73 


M x + 1) 15(r - 1} \ 3tt(_r - 4) ' 70(x -6) 

/(6) = 5.66, maxiourm. at x = 4.5 


A A A 


i 

.1 

vi> 




& 


£ 

fcc 

KK 


https ://t.me/upsc_pdf 


https://upscpdf. com 


https: //t.me/upsc_p(|f 












































Join Telegram for More Update : - https://t.me/upsc_pdf 




o 


^V7'l3: DlVIDED,DiFFERENCES 

The Lapbhge’s formula has the drawback that if another interpolation value were 
inserted then the interpolation coefficients are required to be recalculated. This labour of 
recomputing the interpolaHon coefficients is saved by using Nekton’s general interpolation 
formula which employs what'are called ‘divided differences’. Before deriving this lonHiila, 
we shall first define these differences. 

If (* 0> y 0 )> (x lT y p, (x 2 , y 2 ),../be given pojnts, then the first divided difference for the 

arguments x 0 , x 1 is defined by the relation [x^xji = ^ 7 —^ 7 *.- 

*1 -‘Xq... _ . ; ■. 

Similarly [x i 2 ] = ——— (*2» X J ~ ~ ~ etc ' . 

X 2 -X 1 , * * 3 -X2 

' ' c V r ■ i 1 [xi,x^l - [x 0 ,:g r ] 

The second divided difference foitr 0 , x l? x 2 is denned as i x Q , x 2 — x 0 

The third divided difference forx 0 > Xj,x 3 ,Xj is defined as - .. «■ 


u x x x ] = and sor'on. 

u 0 ->i-^.^ ; ; x 3 -io . ; 


31- 


«s 

g 


3S 

5& 


g 

1 


& 




'-5 

-Vf 

I 


3. 


m 


•u 




>ivv > Ac_cl , cl 01 O^ ) I cx jjU 


^0 
^ i 

^2 


<-f-i r<St 


K ^ 3 

Cm . 


££OCwfI c\.~A 


,«xJJ 

|> < 'M., ] 


ilkgcf cf-'d 


,\,Aj 


iL 


®E£SESSffi?3s:.-7C.i;—- ■ 
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If^PEfJPOtATlON 


Obs-1. TJiediuided differences are symmetrical in their arguments ie. independent of-the. Qrder 


of She arguments ; For it is easy to write tr 0 „ Xj] = 


^o- 


x o *i 


x r“ *6 ■■ 


yo 


y\ 


yz ■ 


(xq - X,) (x 0 - JCg) (x x - roHxx - xf) Crj -x 0 )(xj-x 2 ) 

~ frj, x z , x 0 \ °r [x 2 , x 0 , x i r.and:so Qn. 

Obs. 2. The n/7i divided differences of a polynomial-of ike nth degree are constant. 

Let the arguments be equally spaced so thatXj — x Q ~ x^ Xj = ... = x r —x^j.= k. Then 

- r: lx . x]= yizyo = *y°- ' ' > 

lx 0' xp - ^ ^ £ 


*l-*0 


l*a> 




[xu^l-ixp -Xj] _ X f Ayi _ Ay 0 | 

-X2 — Xq 2h | h h ] 


= an <* in general, (x^ x „- 2 .xj ^ 0 : 

If the tabulated function is a nth degree polynomial, then A’^ will be cmrstaotlHehce the n th 
divided differences will also be constant. ./ - . * 

7 . 14 . NEYmWS DIVIDED DIFFERENCE FORMULA , 

.. y n be the values ofy == fix) .corresponding to the arguments x Qf x 3 ,.. 

x. Then firomthe definition of divided differences, we have to 


fr. x 0 ] = 


y-yp 


so that 


Again 


x~x 0 
y=y Q + (x-jc 0 )[x r x Q ) 

(x,x 0 l - \x 0 , Xj]. 
^> x o.^ =. 


..XL) 


which gives [x, x Q l = [x Q , x x ) +■ (x - Xj)(x, x Q , x^ 

Substituting this value of {x, x Q | in (1), we get 

y^yQ + ix — jColteo'Xj+ix-XQXx-xJiXrXpxJ. --(2) 

, k-Xo.Xjl-Iro Xi,X 2 ] * 

Also lx, x Qt x v x 2 ] =-:- 

X - x 2 

which gives lx, x^ xj = [x 0 , Xj, x 2 ] + (x -x 2 )[x, x£x lr x 2 ) 

Substituting this value of lx, x 0> Xj] in (2), vve obtain 
y =y 0 + & ~ x x)t x o> x il + - x o^ x “ X 1^ X 0< x i> + (f - x 0^ x ~ X^X -x 2 )(x„x 0 , X x , X z j 

Proceeding.in this manner, we get 
y = /Cx) = y 0 + (x - x^fiQ, Xj] + (x - x 0 )(x - X,)fx 0 , X x , x 2 ! 

+ fx - X 0 )(x — XjKx X lr X 2 , x 3 ) + ... 

+ (x ~X 0 )(x — Xj). Xx -«x n )(x, Xq, Xj, - X r ) -X3) 

which is called Newton's general interpolation, formula, with divided differences. 
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NUMERICAL METHODS IN ENGINEERING ^IQ SCIENCE 


Example 7.12. Given the values 
x : 5 7 ' 11 13 17 

f(x) : 150 392 1452 2366 5202 , 

evaluate f(9) r using ( i ) Lagrange's formula 

(££) Newton s divided difference formula. 

SoL UhHere x 0 = 5, x 1 « 7, =11, *3 = 13, 

y 0 = 150,^=392, y ? .=. 1452 , y 3 = 236G,y 4 = 5202. 
Putting* = 9 and substituting the above values in Lagrange’s formula^ we get; 

tb-7)(9-llK9-13)(9-I7) - (9 — 5) (9 — 1I)(9- 17) ; 

k (5-7)(5- 11 ) (5-13X5-17) X (7 - 5) (7 -11)(7-13X717) * 93 


(P.T.U.. B. Tfieh., 2005 ) 



= 17 


and. 


(9-5) (9 —7) (9 — 13) (9 —17) ■ 


50 3136 3S72 

—+-+ —- 

3 15 3 


(11-5) ( 11 - 7) ( 11 - 13X11- 17) 
(9 - 5)<9 - 7) (9-11) (9 -17)\ 
(13-5) (13 - 7) (13 - 11X13-17) 
(9-5)(9-~0 (9 -11)(9-13) 

+ (17-5)(17-7) (17-11X17-13) 

2366 578 

-r + ;r = 810 


X 1452 


*206 


x 5202 


(££) The divided differences table is 


X 

y 

divided 

2nd divided 

3rd divided 



differences 

differences 

differences 

5 

150 






392 - 150 =121 
7-5. 



7 

392 

1452-392 _ 

- n =265 

11-7 * 

265 ~ 121 .=24 
11-5 i ■ 

r-A 

\\ 

3 « 

cl) S2 

n 

11 

1452 


457 ~ 265 =32 
13-7 -■ 




2366 —1452 „ rrT 

-—- =457 

13-11 


42-32 

17-7 ~ 1 . 

13 

2366 


-7b9-4S7 =42 

17-11 




5202-2366 7fl;> 

17-13 



17- 

5202 





Taking x = 9 in the Newton’s divided difference formula, we obtain 

* /(9) = 150 + (9 — 5) x 121V (9/- 5)(9 — 7) x 24 + (9-5)(9-7X9- 1-1) x 1 - 
= 150 + 484 + 192 — 16 = 810. 
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if ded 
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Vt = 1 
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as 


8 

C» 


V? )X ] 

vV 

sa 


rarnTT i "TfrrTB*i 


interpolation 


: .. < 




>/>r - . 


'Example 7.13. Newton's divided differences formula; evaludtcf(8)and f(15) 


£7 pen : 


x: 4 5 7 10 11 13 

f(x): 48 100 .-£94 900 1210 2028 

Sol. The divided differences table is 


(V. TM.M&tfe2008) 


X 

fix) . 

1st.divided 
differences 

- 2nd divided 
differences 

3rd dioided 
differences 

4f)i divided 
differences- . 

4 V 

ff 48 


■ : ' i 





52 



:>v 

5 

100 


15 



0 

. V 

97 


1 . 


7 

294 

202 

21 

1 

0 

10 

900 


2T 


0 



/ 310 


1 


n 

121P 

/ 409 

33. 


* 

13 

2028. 




\ 


~~ 




\ 


Taking x = 8 in the Newton's divided-difference formula/we obtain - - 

fl8) = 48 + (8 - 4) 52 + (8 - 4) (8 - 5) 15 +-(8 -.4) (8 - 5/(8 - 7) 1 
= 448. 

Similarly/(15) = 3150. - 

fecample 7.14. Determine fix) as a polynomial in x for the following.data : 

V-' x;• - 4 -l 0. 2 5 

fix): 1245 33 5 9 1335 . 

(V.T.U., B. Tech., 2007) 

Sol. Th» divided differences table is 


Applyfng Newton's divided difference formula- 

* 1_ A*) = /&o> Oc .+. ... 

_ . - = 1245 + (i +["4y(i-'464) + (x-+ 4) (x + 1) (9,4) 

• +ix + 4)W +■ lXx-0)(-i$ + Cx -i- 4)(i + lMx-2)(3> 

= 3x‘ 1 .—5r?> 6bc 4 ^ ! r4x 


fc .3 

8 

EE 


x 

' fix) 

1st divided 

2nd divided 


4dx divided 

| • 
! 



differences 

differences 

differences 

differences 


- 4 

-« 1245 





r- I 



~ 404 i 




6 [’ 

- 1 

33 


94 





■'* 

- 28 

- 

- 14 


y, |- 

0 

5 


10 


3 

x 



2 


13 



2 

9 


. 88 






442 




t 

.5 * 

\ 1335 




V 

fer i. 
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Nu-waical Analysis 

=.^f(i;6000QO+.50Q600) +2(0.961538 + 0.86206.9 + 0.735294 + 0.609756)] 


0.2 r 


= 0.7839314, 

l ~ 0.78373, correct to five significant figures. 
- ‘ - The exact value 


- [-j*fx = ftan 1 x] 

„l+x 2 l J( 


- tan 1 3 — tan ! 0 


= - = (*7853981 
4 • ' “ '■ 


• - - ~-dx = 0.78540 , 

_ 0 > + * 

correct to five significant figures. 

. Th e error is = 0.78540- 0.78373 

= 0.00167 
Absolute error = 0.00167. 


Example^ Find the value of taking 8 subintervai 

decimal places by Trapezoidal rule; 

Solution.: Here " 

- /M = log 10 x, 

« - l, b = 5, n = 8; 

. f h-a 5-T 

/? = —= ~-= 0.5. 

» 8 


s correct to four 


r 1.0 


1.5 


2.0 


. ---_ 3 : 0 .- _ 3 . 5 . 4 0 ■ 4 -^ 

fix) 0.00000^0,17609 0.301 Oj 0^794 .0.4771 2 -0^407 0a 

- ■ - y » y 2 y. 


510 




Using Trapezoidal rule we can'•tirfite 
h 


.60206'0765321.0.69B97 

^y 7 y s 


7 2 f(7o + y 8 )+ 2 ( 7 , +y 2 +y 3 + y< + y,+ y< . +y 7 jj 

05 f , 

~ —[(0-00000 + 0.69897) + 2(0.17609 + 0.30103 + 0.39794)] 


-■£§j 






is 

s 


. . .__ 


? • ’ 


& 
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Numerical-Integration 


05 ;■ 


[2(0.47712 + 0.54407 + Q.602Q6+ 0.65321)] 


I 

§ 

g 

i 

I 

m 

3 


= 1.7505Q25 
. 5# 

/ = Jiog 10 x/fx = 1.75050 . 

i 


0.6 ' 

Example^V'Find the value \ C: ^ X f taking n ~ 6, correct to fwe significant frg.iires 

by Simpson's one-third rule^_ _ — 

Solution We'have _ V 

: : v f(x) =e x , .. 

a = 0, b ~ 0.6, n =p 6. 

0 . 6-0 _ 01 
n 6 . 

y ' ~Q-Q ~oTi 0.2 ' 0.4 05 ; ~ 0.6 

y =/(x) 1:0000 K105T7 1.22140 1.34986 1.49182.1,64872;$482212 

*o h ' ' *>3 i'-n 


The Simpson's rule is 

/irV 
3 

0.1 


/ 


‘ 


§ff 




= “[(y<> + y&) + 4 (y + y.' -ys) +2 (y2 + yi)] 

= y[(l.0,0000 + +82212) +4(1.10517 + 1.34986+1.64872) + 2(1.22140 +1.4918.2)] 
= ^[(2.82212) + 4(4.10375) +2(2.71322)] 


... 0.8221 186 == 0.82212 : 

= 0.82212. ■ " 

Example 5^Th^. velocity—pf a tram which starts from rest is given by the 
.following table, the time being reckoned in minutes From the start and the 
speed in km/i|b.ur- . : 


i 2 4 

(minutes) 


6 B : T O • 12 - 14 16' I B .20 


16 28.8 40 46.4 51.2 32.0 17.6 8 3.2 0 


(km/hr) - 


Estimate approximately the total distance run in 20 minutes. 


t 4 ; • 
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Solution : 


d$ 

V — — ds = u.df 
di 


f.dt 


Jjs= j"u.< 

• 20 ' ' 
s= ju-df. 

■' '■* 0 
The ,tcain. : star-ts.from; res t y .the. velocity. v = 0 when, f ==./)_. 

.The-given table of velocities can'be written 


/ 

0 

2 

4 

6 

8 

10 ' 

12 


16 

M 

20 

V 

0 

16 

28.8 

40 

46^4 

5l / 2 

32.6 

17.6 

8' ” 

3:2 ‘ 

‘6* ; ' 


Vo- 

: *1 


*3 

. y. t - 

y s 

y* 

y? 

y 8 

- V 9. 

y in 


. -■ -2-V - 1 

A = — hrs = — hrs . 
60 - 30 


The. Simpson's rule is 
• * ' 20 . . 


S = ~^[( y 0 4 ' 2/lo) + 4(y,+y 3 -l-y s ■' y? I y 9 }+ 2 (y 2 +y 4 4- y 6 +y 8 }] 

O ' • 

.i 

-—[(o + 0) + 4(16 + 40 + 515 +17.6 + 32} + 2(28.8.+ 46.4 + 32.0 + 8)] 


30 x 3 


-—(0 + 4 x 128+2 x 1152} 

on ^ 


90 

= 8.25 km. 

/. The distance run by the train in 20.minutes = 8-25 kms. 

Example 6 : A tank in discharging water through an orifice at'.a depth of x. 
meter below the surface of -the^v^ater. whose area is Am 2 . The following are the 
values ojf .x for the^-corresponding values of 


-/l 

\ 1.257 

U39 

: 1-52 

1.66:; 

1,809, 

• V.9.62- r 2,123 

2.295 ■2.462. 

2 ; :$5d;.. -g^:2r7: . — 

x" 

1:50 

.1,65 

.1.80 

1-^5 

.2.10 

2.25. 2.40 

‘2.55 . 2.70. 

.2:85 - a-OjO 


3.0 


1.5 


Using the formula (0:018)T = . f -j=dx ,calculate T. the time in-seconds for the 

J Jx 


level of the water to drop from 3.0 m to 15 m above the,orifice. 
Solution : We have h = 0.15, 


J§ 


ii 

H l 
" 4 %. 


% 

I 

% 


% 

% 

*2 


-- dcc: stsssssec; 




fe 


50 
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if 

• -! 

-'■> a 

• v "Ss 

1 


§■ ® 

1 ! 

*4 


■*£ 


V 


Numental-Ini^aVon 

, the table of values of x and the corresponding values of ^ is 


2.40 »2.55 2.70 2.S5 ' 3:00 


y'= 


71.50 1 . 65 “ 1780 ... 1 . 95 _ 21_°_? ? - 5 

: - "37 ' ' • - 

A 1.081 * 1.132 1.182 1.249 1.308 1.375 1-438 1.498 1.571 1.632 


r. 


1-025 


if. ill 

i l m 


■dfsing Simpson's rule, we get 


fjg 

a tiW. 

:§ 

s jae 

:: c- T-'js •; ■ 

-• O -*§ 5 :V- : 

31 


=^[( y „+.y 10 )+4(3/1 + y.i + ys + y ? 1 y-j )* 2 ^ 1 y< + y<> + y«)1 

■ _ Hl|j[T.ois'+ 1.632.) 4- 4(}-081 +1.182+.1.308 + 1.438+1.571)] 


■ \^2{U32 + 1.249+1-375+1-498)] 


3* 

6 =-19743 


r.s 

Using the formula 


= 1 . 9743 '. 


- •/> 


crafty 



6 


vr 


/ 

vve get i, '/.■?' o ' 

, -f — = 1 tT) sec (approximately) 

■^.018 

/. T = 110 sec.. . 


Ex,mpl«,7 = Evaluate HfV,„ ordi„«. 

.Solution : We ; fiave 
The points of division are 


n .1= 7 - > rr>6e - . . . ; - 


k 


1 2 3 4 5 

o. ~ ~ > 1 - 

6 6 6 6 6 
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By Weddle's rule we have 

Jtr 

2 


N m?.cjical Analyst 


l = 


~ 0.162sin 2 0 d<p 


= 3hr f 

-fo^ + y>2My J+ys+y7+y]] )j 

- 3fa, 

. + I0^ 2+y \ + y s + ' y io) + 6(y 3 +y 9 ) + 2y s J 
3/r 


- 3/r r/ 

" 2 3 ^ (0 " 4559 ' 0 979541 + 0.937283 + Q -. 92 B 322 )] ' ' ' " 

3* r ■ ' * ' ' .... * 

: 6(0-988067 + 0.92829.1) + 2(0.95864S):J ” 

1 --‘1.505103504 . 

S"r„: v H r /w -• --!■...... = 5 , , ati „ s . 

5 6 <<' multiple of 



Jog'Xdx l0 [vo+5 yi +y 2 + 6j/3 + ^ +Sy s +y 6 .] 


= 3.x (02) . 

T"^ + 7 -^^ 1 «« JS 66 +1 ^ +wwd+ 

= 0^#Dl4643j 
— 1.1327858 


5.2 

pog, r=L827858 . 
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Numerical Integration . 

Exercise 9.1 


1. Evaluate by Trapezoidal rule. 


o 




_„ 


2-. E valua te' taking 10 intervals by- Trapezoidal rule: 

o . - 

3. Given that c° =1, e\ = 2.72, e 2 = 7.39, e 3 = 20.09, e 4 -= 54.60^ ^^ 

appr b xirrlatio n : value of j*e X dx' by Trapezo id a I >tdle. ' ; W 

~~ - o / 

i 


an 


ri u 


4. Evaluate Jvl -* 3 dx by (i) Simpson's rule and (ii) Trapezoidal rule, talcing 
six interval correct to two decicriaI places. 

i _^ 

5. Evaluate jVsinx dx taking x = 6, correct to four significant figures by 
(i). Simpson's one—third rule and (ii) Trapezoidal rule. . 

. . Ml V f . 

6^ Evaluate J talcing 4 sub intervals, correct to five decimal places 
'<( - ] x . 

(i) Simpson's one third rule (ii) Trapezoidal rule. 

_ i 

7. Compute by Simpson's one-third rule, the integral |x 2 ;(l —x )dx cdrf^cj- fo 

• - -- ..■ o c ^ 

three places of decimal, taklng step length equal to 0.1. 


0 ’ 


m 


i 


B. Evaluate Jsinx 2 dx by (i). Trapezoidal rule and (ii) Simpson's one-third 


rule, correct to four-decimals taking x ~ 10. 

3 


9. Calculate approximate value of jsinr dx by using (ITTrapezoidal rule 

-3 

and (ii) Simpson's rule, taking n = 6. 
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10. Find the value of J VcosxJx by (ii) Trapezoidal rule and (ii) Simpson's 

o ... 

one-third rule taking x = 6. 

1.5 

\r Compute .^ X dx by (i) Trapezoidal jjuleand {.ii) Simpsontsuone—thirdrule 

taking x =10. . 

0.5 




£ • 


12- Evaluate [- rdx . taking ii =10, .1^ 6) Trapezoidal rule and 

/ .0 * 


(ii) Simpson's one—third rule. 

0.4 , . ' ; .vb . 

13. Evaluate Jcosx<ix taking four^equal intervals by (i) Trapezoidal rule and 
(ii) Simpson's one— third rule. 


A 



1-4. Evaluate JVcoS.vi 

o 



0 

w 

A 



1 2 . o 
rx -f- Z 

15 Evaluate 0 

/ x 2 4- 1 

3 



0 

n = 12. 

N 

1 

:-.X 


2 i 

N 



16. Evaluate f- 

© © 

-° r *r 


n 1 + x 

CO 



V 

Gh 

m -X4 

: V: 

1.6 

17T Evaluate f*~- 

J slnhx 




— 0.4 


1 


five decimals. 

-* ’ / * 


dx taking thirteen—ordinates by Weddle's rule correct to 


* i _- 

18. Using Simpson's, rule evaluate sin jmTx with seven'ordinates. 

. ’ 0 
2 - 

19. Using Simpson's rule evaluate J^x —1 Yxdx with five ordinates. 


$ 


a 

S 

:f 

.3j ; 

3 I 




g 

c: 


6? I. 




s 

tXs 

r,4 

g 

i 

& 

i 


£2 

1 

1 






sc* 

Ic 

i 
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Ord inary Differential Equations 


Example 8.1 Using Taylor's series method, find the solution of the initial 
value problem 

— dy 

-jy = t+ y> xi)=o 

< / 

at / = 1.2*with h = CkL^and compare Ahfi result with the closed form solution- 

Solution Uet us compute the first "few derivatives from the given 
differentiafcjequation as follows: 

y’=t+ y, y" = 1 + /, /" = y", / v = y"\ . y v = y lv (l) 

Prescribing the initial condition, that is, at t 0 - 1, y Q ~ y(/ G ) = 0.,-we have 

yd = 1 

Now, using Taylor's series method, we have / 


ya = 2. j'o" = > 0 V = Jo ~ 2 


jKO - yo + - l o)y'o + ~—~~ L ’Jo * ~— ~~—y'o" 

L o 

. (* ~ ^o) 4 IV (f — ^o) 5 V . 

+ - y 0 + -— y 0 + — 

24 70 120 


( 2 ) 


Substituting the above values of the derivatives, and the initial condition, into 
(2), we obtain 

Q01 0-001 0.0001 , 0.00001 ■ 

y(l.l) = 0 + (O.D(l) + -±—(2) + ~ (2) + —r(2) — C) + --- 

2 6 ^24 120 

£ 

■ 0.001 o.ooo^;- /;0,(S)001 

= 0.1 ,-O.Ul + T +_ ^^p + - : 

= 0.1 + 0.01 + 0.000333^46%0p^§ + 6.0000061 + — 

= 0.11034 14. 

Therefore, 

=ithl(J34 


03414 = 0.1103 


v C1.1.) = /, , 

•4? 

Taking y, = 01 103 at r = 1.1, the*values of the derivatives as computed from 
Eq. (1) are 

y{ — 1-1 + 0.1103 =1.2103 
y," = 1 + L2103 =2.2103 

y\~ = yf v = >7 = 2.2103 

Substituting the value of v, and its derivatives into Taylor's series expansion (2) 
we get, af(ef"retaining terms up to fifth derivative only 

v , ’ 0 ~ fi) „ 0 - fi)‘ ,„ t (t ~ h) iv 0 — fi) 5 v 

2 6 24 120. 

= 0.1103 + 0.12103 + 0.0110315 +0.0003683 + 0.000184 + 0.0000003 

= 0.242934 I 


Cf>& 


1 


t 

T 

4t 
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Numerical Methods for Scientists and Engineers 


Hence, 

y(\2) = 0.2429341 = 0.2429 

To obtain the closed form solution, we rewrite the given IVP as 
dy 


dt 


-y - x or d(ye ~•) = te 


On integration, we get 

y - -d (fe~ l + e~ { ) 4 - ce f — cd ~ t— 1 
Using iheLinitial condition, we get 


0 = ce — 2 or c = ~~ 
e 

Therefore, the closed form solution is : .g, ? 

when t - 1.2, the closed form solution becomes \ 

>(1.2) = —i.2 ■- 1 + 2(1.2214028) - -2.2 + 24428056 = 0.2428 


( 3 ) 


( 4 ) 


Comparing the results (3) and (4), obtained numerically and in closed form, we 
observe that they agree up to three decimals. * 


e 
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Computer-Oriented Numerical Methods 


EXAMPLE 7.2 Given ^ = 3^+- and y(0) = 1- Find the values of y(0.1) 

dx 2 

and y(0-2), using the Taylor senes method. 


Solution 

Given' 


y' = 3x + yi2 and x Q ~ 0, y 0 =■ 1 


Taylor series is given by 


=x 0 + h yo + —y'o + —yo+~y'o +-- 


3! 


(i) 


To.find y(O.l): 


, V v , o ■ -Vo 

y =3*+- => y 0 =3x 0 '+—■ 

■.3(0) + i 

= 0.5 


= 3 +- 


Vo = 3 + 




= 3-^ 
2 

- 3.25 


y = y -> 


y =— 
2 


>o 


2 

3.25 

2 

1.625 

//' 

Zo. 

2 

1.625 

2 . 

= 0.8125 




Equation (i) => 

y | =I + ^0.1)(0.5)+Y_ (3 - 25) + 


(0.1) 


(1.625) +^L (0.8125)+.... 
6 - - - 24 . • 


- J.0665 

yi =><0.1) 

= 1.0665 
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To find y(0.2) : 

Taylors Series formula for y( 0 . 2 ) ; s 


Dm*** 


Equator 


& -y\ +hy' +—y" + h 4 iv 

2! 3i >l f i+•- 


^~ 3jr+ 2 ^ yl=^+~ 

z . 2 

= 3 ( 0 . 1 ) + 


i-0665 


: 0.83325 


00 




>'~ = 3 + - 


O 



Equation (ii) : 


>r = 3 + i 
2 


= 3 + 


0-83325 


- 3.416625 


y 


=* >r=^- 


y^ 

2 




3.416625 
~~ ' 

= 1-7083125 

-yr 


_ 1 . 708 , 3125 
~~2 ~ 
= 0.85415625 


* 1.0665 + (0.1) (0.83325) + (3.4 16625) 

, (0-1) 


6 


= I- 167196 . 
?2= X0.2) 

= 1 . 167196 . 


EXAMPLE 7J Ob,*, ,(4.2, ,„ d + 4) 


given 


dy_ J 

dx r 2 + „ ’ y(4) = 4 - ' akjn S h = 0.2 
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Computet-Oriented Numerical Methods 


(z’X. 


Solution 
Given y r = 


x 2 + y 


. Xo = 4 ^0 = 4, x, = 4.2, x 2 = 4.4, A = 0.2. 


Jaylor series is given by 

To find y(4.2): 


yi=yo + hyo + yo" + 3 , fo" + •* 


(0 


y = -“2 - 

x + y 


y 0 ^- 

* + y<> 


l 


4 2 + 4 


= 0.05 

y" = -1U 2 + y)“ 2 y r => >o' 
. - 0.05 


">o 




( 20 ; 

= - 0.000125 _ 

' = 2{x 2 + v)- 3 (y') 2 + (-& (x 2 + y)“ 2 y" 


2 W>‘ few 


(•*0 +• 


2(0.05^^J$#5 

= -- 


(20ra5f«£ffi0r 


= o.ooooot#3?£ 


' a 


/. Equation (i) => 

y, = 4 + (0.2) (0.05) + 




^' 2 - L (- 0 000125) + £0.0600009375) + - - - 

2 6 


— 4 + 0.01 - 0.0000025 + 0.00000000125 
= 4.0099975 


y,~y(. 4-2) 

= 4.009998 


To find y(4 A): 

Taylor series is given by 


i , h * h 

y-i = y\ + H + 2 ) )’i + 37 + 


(ii) 
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Numerical Solution of Ordinary Differential Equations 


y = 


x 2 + y 


> i ="2 - 

*! +Vl 


1 


(4.2) 2 + 4.009998 


71 = 


= 0.046 ia9 

r 

y\ 


-0.046189 


(4 + y, ) 2 [(4.2) 2 + 4.009998] 2 

- 0.000098542 


yr- 2W> " 


^1 


tf+y,) 3 (*? + y >) 7 


2(0.046189) z ^ + ._ 0.000098542 


[(4.2/ + 4.009998] ((4.2/ + 4.009998/ 

= 0.000000420469 + 0.00000021024 
= 0.000000630904 


/: Equation (ii) : 


y 2 = 4.009998 + (0.2) (0.046189) + 


(0.2) 2 


(- 0.000098542) 


(0.2) 3 


(0.000000630704) 


= 4.009998 + 0.0092378 - 0.00000197084 + 0.00000000094094 
= 4.019234 


EX A MPLK 1A Find >(0.1), given y' = .\ 2 y — 1 , y( 0) • 1 . 
Solution 

Given y*-^x 2 y - 1, x 0 = 0, y 0 = 1 > .x l = 0.1, h = 0.1 
Taylor series is given by 

h 2 /i 3 

= >o + h y'o + — y'o +— yo" + -•• 
y'=^y-i => yo=x 0 2 y 0 -i 
=10(1)-l 
= -1 

f~ = xry' + 2xy => 7<f--'■o 2 >o' + 2x 0 >' 0 
. =0 


(0 


y'" - x~(y') 2 + 2xy' + 2xy' + 2v 






m 

m 




fc. 
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Computer-Oriented Numerical Methods 


y'v" = *o (yi) + 2x 0 y' 0 +2x 0 yQ.+ 2y Q 
= 2 


/. Equation (i) = 


> i =l + (0.1)(-i) + M 7 (0) + M. (2) + ... 

Zl ■ . 6 


. = 1-0.1 + 0.0003333 + ... 

= 0.900333 
y,=y(0.1) . 

= 0.900333 


7 . 1.1 


Differential Equations 

The simultaneous first-order differential equations of the form 


dx 


and 


dz 


dx 


= / 2 (.w) 


with initial values y(.v 0 ) = >'o and z(x 0 ) = z 0 

To solve this system of equations at an interval h , the increments in y and z axe 
obtained by using the formulae ; 


- ; , * ^ w 
,vi = v 0 + *y 0 + x[ y o + 37 yn + - 


and 


h 3 


Z| =_z 0 + ^'Zo + —Zo+ —Zo" + -- 


EXAMPLE 7.5 Solve the differential equations using the Taylor series • 


- ■ : dx dx, 

given that x — 0, y = CL z == 1 1 - 
Solution 

dy 
dx 
dz 


dy dz y 

— I + xz, — = - jcy, for0/= 0.3 


; = l + JCZ 


dx 


■-xy 


x = 0, y == 0 ? z = 1, /z = 0.3 


m: 




£ 

£ 


I':':':- 


1 


li 


I 


L- ■ :: 
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Numerical 'Solution of Ordinary Differential JEe, 
Taylors series for y' is given by, 

h 2 h 3 

y { =y 0 +A %+~>'o +~ y'o +•- 

y'=l+xz => yo' = i+x 0 z 0 

= i+(0)(l) 

= 1 

. .y" = xz' + z => y 0 "=x o zo' + zo 

= V-^yo)+zo 
- 0+1 
= I 

/" = x(zy 7 I z' + z' => Jo" = a 0 (zo'f + 2J 
■ =x 0 (-jc 0 y 0 ) 2 + 2(x 0 y 0 ) 

= 0 




3£ 

! I 


(i) 


Equation (i) - 


( 0 . 3) 2 


0) 


(0.3) 3 


( 0 )+.:. 


y? = 0 + (0.3) (1) + 

= 0.3 + 0.045 
= 0.345 
V) = >'(0.3) 

= 0.345 

Taylor series for z ' is given by 

11 j t 

z, = + hto +-Tz'o + ~y- z"’ + ... 

z'--xy =* Zo=- Xoy 0 = 0 
z" = -(xy'+y) => = -Xoy 0 ' + y 0 

= -(0)(l) + 0 
= 0 

_ ^ y + y + y} =$ Zq = ■ _ -*b3'b 

( 0 ) ( 1 ) + 2 ( 1 ) 

’ = 2 . 

Equation (ii) ==> . ~. r 

z, =l+(0.3)(0) + ^-(0) + ^/<2)+... 
2 6 


oi) 


I 

& 


f; 








v.<: 


■ - 
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= 1 + 0.009 
= 1.009 
z^zCO.3) 

. = 1.009 

EXAMBLE 7.6 Find >(0.1), >(0.2), z(0.1), z(0.2), given 

•SZ = x + z , — =x-y 2 and >(0) = 2, z(0) = 1 

Given 

y ' = x + z, Z' = X - x 0 = 0, y 0 = 2, Zo “ 1, h = 0.1 

Tofindy(O.l): 

"Taylor series is given by 

.2!-' v 3! 

y'. = x + z, z ' — x-y 1 

y" = 'l + z\ z" = I -2yy' 

y'" = z" 

t" — — (2>>~+ 2(y') 2 ] 
y' = x + z => y 0 ' = x 0 + z 0 
= 0+1 
= 1 


>i = yo + ty> + —70 +— 7 o +- 


© 


y" — 1 + z' = l + x - y 2 
= 1 -2 2 


yo" = i + ■ yo 


Hquation (i) : 


y 1 = 2 + (0.1)(l) + 


(0-1) 2 


(-3) + 


(O.I) 3 


2 ' ■ 6. 

"=2 +0.1.-0.015-0.0005... 

= 2.0845 

Tofindz(0.1): 

Taylor series is. given by 

• /i 2 „ A 3 ~ 

z, = z 0 + Az 0 + — Zb t 

zo'^o-yo 2 
= 0 - 2 2 
= -9 


(-3) + 


(ii) 
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Numerical Solution of Ordinary Differential 

= 1 ~ 2y 0 
= 1 - 2(2) (1) 

= *-3 

2T'=-P^+2(y5) 2 ] 

~ — 12(2) (- 3) + 2(X) 2 J 
= -f-I2 + 2] 

= 10 


.*. Equation (ii): 


Z, =1 + (0.1) (— 4) + 


( 0 - 1 )" 


(-3)' + 


(0.1) 3 


2 6 
= r 0.4 + (- 0.015) + 0.0016667 
= 4 - 0v4 - 0.015® 0.001667 
= 0.5-86667 


(10) + .-.. 


To find 

Taylor series is given by 


. , A 2 „ IT 
y 2 = y, + Ay, + — y, + — y, 
2 ! .>! 

y{ = *\+z\ 

= 0.1 + 0.586667 
= 0.686667 
y\ = > +zi' 

= 1 + X) - 7i 2 
■'= 1 +0.1 -2.0845 2 
' =-3.24514 

M // 

3»i = Zj 

= 1 - 2y,y,' 

= 1 - 2(2.0845X0.686667) 
•=- 1.862415 


(iff) 


Equation (Si).: 


y>= 2.{S®5 + (O.l) (0.686667) +' 
+ ^^-(-1.8627)5) + ... 


( O .)) 2 


(-3.24514) 




g 

S3 


SC 

e 




© 
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Computer - Oriented 'Numerical Methods 


= 2.0845 + 0.0686667 - 0.01623— 0.00031045 
= 2.13663 

To find z(0.2): 

Taylor series is given by 
' , A 2 ,, h 3 

z 2 =z l +hz' l + — z: + ~z? +... 

Zi =x l -y i 2 

= 0.1- 2.0845 2 
= - 4.24514 

Zi" = 1 - 2 y x yi . 

_ 1-2(2.0845) (0.686667) 

= -1.862715 

z[" : = — [2>, yf+2(?,') 2 ] 

• =-[2(2 .0845) (^3.24514) +2(0.686664) 2 ] 
= -(- 10:943023] 

= l2.585^i^?<o . 

.-.Equations vVs^’fo^ 


^ c \ -r o 
/ 'Z v---- ^jg» 


Z2 = (0./r+((Q5^4.24514) + ^L(- 1.862715) 

2 ; 

+ % (12.585967) + ... 

= 0.1- 0.424514 - 0.009313575 + 0.00209766 
= -0.33173 

7.1.2 Taylor Series SViBtnod for SoG^nd-orcSsr DitfsrsritiaS 
Equations 

The differential equation of the second-order can be solved by reducing it to a 
tower-order differential equation. A second-order differential equation can be 
reduced to a first-order differential equation by transformation y' = Z- 


' $ 4 T 


Suppose 

Le_. y"=M'y,y') -= ■ 

is the differential equation togetherVwith initial conditions. 

y(* o) = yo and /(*o) = yd 


(7.2) 


(7.3) 


( i 

! i 


£ 

u. 

k 


y 

fc 
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_.r\ 


-ll '. i -4 


where y 0 , y 0 ' are known values. 

Setting /> = y', we get y^ = p-' 

Tire equation (7.2) becomes 

/? y 7 p) with initial conditions 

> W = J?q 

/«=y 0 ' 

where )%-y{/ are known values. 

B> patting y ' =7>, y" = p\ Eq. (7.2)becomes 

p'-My^p) 


CM&- 

P-5> 


with initial conditions 


yM=yo ''i 

.V/C*©) = V => p(^-o) = Po 
Solving p ' by using Eqs. (7.4) and (7.5), we get. 


j2 j -3 

■ft = ft + APO +—Po+^7 Po' ■+ 


P.6X 


where /z = x } - x 0 

Since p ~ y\ we get Eq. (7.6) as 


^2 ,3 

* “ft +/ ^+— yo+~yo+- 


Stmilarly, proceeding in similar manner,, we get 


/,* 


>2 = y. + ^ +— y" + —■ yr +. 


We calculate y.,, y 2 • - - 


EXAMPLE 7.7 Solve y" = y -t- given yfO) =;1, y'(0) = 0 and calSSate-^Otj), 

Sfohttion „ £ 

• - ■ < . * 

Jo' = 0" • • 

y^j+jy' j. 

^ ~~ yo" ~j$+xiffo 

-1t(0)(0) 

v =r".7 ■’ -■ * 

differentiating with respect to x 


/" = y'+y'* xy~= 2/ + xy' 






rrvrrTrm-~-ni ! 






a ss SR^/:-a-.gK^y«g L a; 




p&s 




I!* 

a 


1 


K£ 

S5 

14 


sc: 

r 




p 

1 
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Computer-Oriented Numerical Methods 


H 


yo" = 2 y't> +x o- v o 

= 2(0) + 0(1) = Q 


y iv = 2y" + y" + xy"' 


~3y" + xy'" 
y‘ v = + x 0 y 0 

= 3(1) + (0) (0) 

= 3 

y* = 4y'"+xf 

/=4 yv+xv-yo” 


— 4(0) + 0(3) 

=o 

We know that Taylor series is given by 


h\ 


y(x) = yi = >o + h y'o +■^7% + -li y ° + ai + "* 


3! 


4! 


(0.1)* m . (0.1) 3 _ . (0.1)- 


= i + (0. i)(0) + (i) + , 

2 . b 


TO) +; 


24 


(3) + -- 


= i+ ,(ai)l + (o^ ( 3 ) + ... 

2 24 


= 1.0050125 

EXAMPLE 7.8 Find y(0.2), given y" + y = 0, >(0) = 1. y'(0) = O. 


Solution 

y" = - v. xq = 0. y 0 = 1. y 0 ' = 0,h = 0.2. 

To [mdy (0 2) 

We know that Taylor series is given by 


}l n 

y ( = y 0 + h y'o + wr y'o + ^T y o + - - - 


<*) 


■y =-y 
/// - * 
y = -y. 


y iv = -y' 


y = - y 


’=> yo"=-yo = - 1 
=? * y 0 "'=-yo' = ° 

=»- yo h =- >'o"=-(-!)= i 
■ => yo v = -yo'"= 0 


Equation (i) 


y, =1 + (0-2)(0) + 


(0.2) 2 


(- 1 ) + 


(0.2) 3 


( 6 ) + 


( 0 . 2 ) 


24 


0 ) + — 




24 
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= 1 -fc© + &(MM667 

= 0L9SO0667 . . ■ - 


K01S 




. ■ ■ .... 

7..1 Solve — = x + >, given y( 1} ~ 0 tod get y( 1.1), y{ 1:2)' &y Xay& s<mes; 
ax 

method 

{Aiis. >(1:1) = 0.110342, >(1.2)■= 024281} 
7.2 Using the Taylor series method, find correct to four decimal places the 
valtre of : >(0.1) v given - ' r - ' 

dy o 2 # 

-~ = * !/ aod>(0) = 1. 
d-r 

(M. >(0.I) = 1.11145}. 

73 Using the Taylor method, compute >(0.2) and >(6:4) correct lo four 
decimal places, given 


— = 1 - 2^> and >(0) = 0. 
eft . 

.[to. >(0.2) = 0.19475, >(0.4) = 0.359384) 

7-4 Given- ■— = 3.x + — andy(Q) = I* Fkid the values of >(0.1) and y(0.2), 
dx 2 

nsing the Taylor series method. 

(/ter. >(0.1) = 1.0665, >(0.2) = LI67196) 
7.5 Solve, by die Taylor series method of third-order the problem 


“ = (x + xy^e x ,#) = ! 
dx 

Co Tmd >, for x = 0.1,0.2,&3. 

[Ans. y(f)J> = 1.0047, ^(0.2) = 1,01812, >(0.3) = 1.039953 
7.6 Solve by the Taylor series method (of fourth-order) 

— — ^ +1, >(0)= 1 at x = 0.2, 0.4. 

— dx 

■ £4/^. >(0.2) = 1.226, >(0.4) = 1.54205} 
7-7 l£stn« the Taylorseries mclhodt solve 

— =\x 2 — y t yQ 0) = 1 al.r =0:1,0.2, 0.3 and 0.4. 
dx 

[Ans. >(0-1) = 0.9052, >(0.2) = 0.821-3, 
>(0.3) = 0.7492, >(0.4) = 0.6897} 




\ 






\ ■ 
l;.. 


<y, 

*S 


% 

% 

% 
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dy 


7.8 Find>>(0.1), given — = x + y,y( 0) = 1. 


7.9 Find >(0.1), >(0.2), >(0.3). given 


[ Ans. >(0.1) = 1.1103) 


, X s +xy 2 , m _, 

y =-:-,>( 0 )= 1 . 


[Ans. >(0.1) = 1.0047, >(0.2) = 1.01812, 
>(0.3) = 1.03995] 


7.10 Solve ^r = y+ xi . forx= 1.1, 1.2, 1.3, given>(1)= 1. 
dx 


[Ans. >(11) = 1-225, >(1.2) = 1.712, 
> >(13) = 1.874] 

7.11 Find >(0.1), >(0.2), z(0.1), z(0.2), given ' .* 


— ! =x + Z,—= jc-> 2 and >(0) = 2, z(0) = 1. 
dx dx 


, .{Ans. 2.0845, 2.1367, 0.5867, 0.1550] 
7.12 Evaluate x(0:1), y(0.1),.t(0.2),>(°- 2 )*8 iven . 


dr y 

given x~ 0 , 7 = l at t — Q, 


dx , dy 

— ~ C, + 1 » — — — l x 
dt 


[Ans, jc(0. I) - 0.105, 7(0.1) = 0.9987, 
*(0.2) = 0.21998, >>(0.2) - 0.9972) 


7.13 Find 7 ( 0 . 3 ), £(0.3), given- 


dz dy 

— = ~xy f — =1 + xz 
dx dx 


where 7 ( 0 ) = 0 , z( 0 ) = 1 . ■ 
7.14 Solve for * and y 


[Ans. 7 ( 0 . 3 ) = 0.3448, £(0.3) = 0991] 


dx dy 

-— -x + y + t, — = 2 jc -r 
dt dt 


giyen x = 0,7 = 1 at / = 1 . 


- 5 3 

Ans.x = 2c + r + -C + 
6 


3 2 2 3 

r +-t + 
7 2 3 


I 
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of Q'rt&uzry Differential Equations 


7-35 Solve jr«inctj®aE^ ihsTayloE senes method fmd approxbaatc., 
values of ^ r aiid^cQ5jcspamfeg foy.yr. 0._2 given that 


-ym= Vz£ 0 }~ I wA~- — x+.z,~-~x~ y 2 

- dx dx 


E 4 m, y(GJ>= 2 . 6845 ,z(O.l) = 0.5867 
= 2 . 1367 , z( 0 . 2 ) = 0 . 15497 ): 




136 Find t&*Y£dae andy(L2} from ^~ + y 2 —= x 3 , y(l) = I, 

: £&* dx 


y *{ 1 ) ~ ] fay imn&.tB&jT&ylm series method. 

. fXhMftl. J) =. n^^l. 2 ) = 1 . 2015 } 

v 2 ‘ 


XI7 Given — jiry^ =Tfc-withy (0^ = 1 ,y '{ 0) ~0,obtain the values 

■' dx 2 / \d%J - ' A .... 


«o£y(&T) airdy(&2}£correct to 5 decima 1 places, using the Taylor series 

; 

method. 

. H>iy’ y( 0 - 1 ) = 0 . 995 , y( 0 . 2 ) = 0 . 9 . 81 ] 

X3S Using the Taylor.series method, find y(&. 1 ), >(0.2), given y" + xy ~ 0 
- aod y(G) = % y ^0} = Ol-S. 

{A-ns. >.( 0 . i ); =, 1 . 049 a, v( 0 . 2 ) = 1 . 0986 ] 




7 : 4 s$ 
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Numerical A mtysrs 


ample 


dy 


: Solve the equation ^ = \ -y, with the initial conditron 


* = 0X0.2, 0.3^ EUlGr ' S alfe ° ri,hm tabulate the solutions at 


Solution : Given ^-'-l-y, . with The initial condition i = 0, y = 0 


we have 


fix, y) ~ 1 — y 


h - 0.1 


x « =■ y 0 = 0 


•v. - .r., + /i = 0 + 0.1 =0.1 


Taking n =0 in 
v/e get 


■ r 2 = 0 2, ; lf ^0.3 


y» i'l - y„ - h H\,, yj 


- =y«. + w* 0 . y„,* 

= 0 + (O.no -0) = 0.1 
y, = 0-1 H^y(O.l) = o.l, 

y 2 =y, + h f(x v y ,; 

y 2 = 0.1 +(0.1)0 -yl) 

= 0.1+ (0.1)0 - 0.1) 

= 0.19 

y 2 =y(0.2) = 0.19, 
y 3 =y 2 + h f(x 2 , y 2 ) 

■r. y, =0.19 + (0.1)0 - y 2 l . 
\ = 0-19 + (0.1)(1 - 0.19) 

= 0.19 + (0.1)(0.81) 

= 0.271 

y 3 =y(0.3) = 0.271. 


, s r' 


* Solut ion by Euler's method 

o “ (T ~ 

oi ,0.1 

°-2 . 0.19 

0.3 0.271 


IT 



Nu jr - 
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fjjifmcncal Solution of Ordinaiy Differential Equations 







— 


. ss 




Example ■: Given — ~x 3 -hy , y.(0} — 1 , compute y(0.2) by Euler's method - 
dr 

taking h = 0-01/ 

| SoluUon : Given : ^ 

«fjC 

with the initial condition yfDJ = : :t .3 

We have - " 

' /(-*, y) = * 3 + y 

,. x 0 = a, y 0 = 1. * = °.°1 

Xj =x 0 .+ h —0 + 0.01 = 0.01, 

x 2 = x 0 + 2/i = 0 + 2(0.01) = 0.02. 

Applying Euler's formula.we get ^ _ - 

y, =y 0 * hfay yjffisf 
y\ = 1 + (O.OD(^ 0 X^ 

= i + (o.ow®f#A. : # 

= i-0i A, 

y, = y(0.01) = 7 #£lf£ 

y 2 =y,'+¥(-^,) 

= 1.01 + <0.01-){* T \t y] _ 

= 1.01 + (O.Ol)KO.Ol)- 1 + 1.0,4] = 1.0201 

.-. y 2 =y^0.02) = 1.0201.' 


Example : Solve by Euler's method the following clifferential equation 

. ... rfy _y-x ..... .. . 

x = 0.1 correct to fou'F* decimal places with the initial coaditioh 

' r ax y + x ' 

y(0) = 1. : 

Solution : Here 

rfy ij-x 


dx y + x 

X/ + X 


the initial condition is 


“I. 


£ 

t£» 

I 

£ 

& 
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ft 




s/ttsse 


= ft.02,- 
= &94, 
2 3 - &&5, 

^ = ft.es, 

*5 = ft-}, 




^+4^11 


>~ J D 


1 ”T > vl.+ 0/ 


= i.O20ft 

= i.&m,. 

— ^r3*t®« =jf. * hf(x v y>> 

\ 


=y,+A|^L-!i 

Vy.+v 


= 10200 4 (0M2)f 

}-ft20\J: 4s-^gggiisl ! 

.. =■ 3-6392 

S'® 8^=1.0392, 









4§L_ y* 

■*" Jtp' ■ 


=S85gggs 

&s= ggpp^ &aa^- 


>j 


/ . 
I 




£ 


?*:* 


i 

g 
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IS 


VJ 


Numerical Solution of Ordinary Differential Equations 
y s =jr( 0 . 1 ) =y 4 + hf(x r 

-y* t h 




\ 

y A~i 


Si' 


a 


= 1 . 0756 +( 0 . 02 ) 


vy* + *<7 

1 . 0756 - 0.08 
1.0756 + 0 . 08 . 


= 1.0928 
y(0.1) = 1.0928. 


- :Vr. - - 

/Example _ : Solve the Euler's modified.method.The.following differential 

dy 2 

equation for.x = 0.02 by taking h == 0.01 — = * V when x = 0. 


Solution : Here we have 


/(x, y) y 

;i = 0 . 01 , X 0 -J), y Q = y( 0 ) = 1 
x, = 0.01, x 2 = 0.02 


we get 


• y, l0) = y 0 + h f (x o' yJ 

= 1 + (O.Ol)(x 0 2 + y 0 ) 

= 1 + (O.Ol)CO 2 + 1 ) = 1.01 

.•.y, (0> = 1.01. 


- 


o -+ * V’ J ' 


Applying-Euier's modified formula we get 

/( x o'yo) + /( z i'y! 0) )j / 
+ 1 +(o.oi) 2 +1.0l| 


w h 


0 . 01 , 7 

= i + — lo 2 


2 

■. 4^,01005 
/. yi (I) = 1 : 01005 , 
h 
2 

0.01 


y, (2) =y 0 ^“ /(^o^o) + /( x vyf ) ) 


j 


=1 +■ 


2 

. = 1.01005 
y c> = 1 : 01005 . 


0 2 + 1 + ( 0 . 01) 2 + 1.1005 




v --— v 
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£5 
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■ 

■ 

& 

vt 

0 


-a 1 


/R-v- 


--3^!’=^® = i-.DiasSi 

, = 1 . 01.695 > ®mkg*-3$ 

0 

4 l)= ?»* f M*>*^ J*4^slSi 


= 1.01065 + 


= 1 . 020294 , 


O-OJR 



-y, ^ 

= i.oi+ 

-:-y, = 1.020204 

^ C. 

~~ y 2 . = y(0-02) = 1 j02tO2£4 






&£= 0.1. 




’»'=*■ 


= -*■ 


series xa&t&crc^ 




^® sa#ea tha ‘ ~=* ***.■&&=% 


SfSSfe — = 3*+ y 2 . 


dy - . .. . 

dx ' '■ y = .?«.»&»•»==. ^..BB31saa£Ea%^l2g£^. t =• 

s^ies method: 


Taylor's algorSStea;•»rf=j?*^ 3 Jgffi ? =ffi'fefe. a4l 8S 
“gg»roxirnaHons to y ($.;»>; C^a»F .^e^aaigsfefiiasis'ap{5 3 rJeciajafe. 


y {l) by Euler's rrM^»<i ^iKa ^dSafecenSal ef •' ^ S 


equation 


<£* ISJftas 


y(0 ' 3) 2 ' €qbvct * a^6b-'fi«a (ftteSBBa^-p&ees takmg step 









IS 




tv 

i-=; 


1 


IS 

te 
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© 

© 

i 
■ ?? 

'* • . dx 

upto two 

tw 

f/i 

$5 * 

9. Given 

fy 


$■5 

If 

y for x 

dx 

= 1 




I 

si 

F i 

I 


Use modified Euler's method with one step to find the value of y. at _t = 0.1 

dy _ 2 ~~ 

to five significant figures, where — x + V , y = 0.94, when .y ~ 0. 

dx 

11. Solve y* — .t ■ ~ y^, by Eulex's method for x ~ 0.2 to 0.6 with /;/= 0.2 initially 
x = 0, y = 1. 

12 r Determine y(0.02), y(0.04) and y(0.06) using Euler's .modified method. 

URSWgFS ““ 

1. 4.0098 2. 1.2375 3. 1.12725 4. 1.052 5. 1.2632 6. 1.01871 7. 1.0202, 

1.0408, 1.0619 -8. 2.42, 2.89 9.1.0928 10.1,039474 11. y{0.2) = 0.8S42, 

y(£>.4) = 0.7(798, y{0.6) = 0.7260 12. y(0.02) =. 1,0202, y(0.04) = 1.0408, 

y(0.06) = UM9 


OV ? 


Numerical Solution of'Ordinary Differential Equations 

6. Find y(4.4), by Euler's modified method taking h = 0.2 from the differential 

dy 2 — y^~ . - - 

equation .— ^ , y ~ l,»wheh x = 4. : 
dx 5x % : 

dy 7 

7. Given —=or+y, with y(0) — 1, evaluate y(0.02), y(0.04) by Euler's 

dx 

method. 

' *r- , , . /n , x ' ■ . , , 


with boundary condition y(0) = 1, find approximately 


set 

MS# 

E 


vSt 
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: ^ se Runge — Kutta method to approximate y when 4 = 0.1, given that 

y = I # when 

rfy 

.r = 0 and — = x + y . 
dx 


f{x, y)—x + y, and h = 0.1: 

’-A jr iryo )=jr o + ^o ==0 + 1 = */ 



ive get 


k i =h /(^ y 0 > 

= 0.1 xi 

= 0 - 1 , 

+ 2 ,y " + ^j 

= (0.1)(/(0 + 0.05), 1 +0.05) 
= (0.D/(0.05,1.05). 

= (0.1 X0.05 + 1.05) 

= 0 . 11 ; 
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+ y) 

= ®-l)(/(O + 0.()5), I + 0.055) 

- 1.0S5) 

=©.i)u;io5) 

# - * • 

=(&1T>G5; 

%—>-&(>+ Ky 0 * 1 3 ) 

=(ai)/(o + o.oi; i + 0 . 1105 ) 

= {ai)/£kl, 3.11.05) 
-(0.1X1.2105) / 

=0112105, ' 

.-:4 f =i(i 1+ 2* 2+ 2i 3 +iJ 


= - (o.l 4 0.22 + 0.2210.+ 0112385$ 
6 

= 0.11034. 




*,=*„ + fc = 0 + 0.1 =0.1 
^ = ^ + 2 ^= 1 + 0.11034 
= 1.11034. 




' zcfeing Runge-KuHa 
«% . ^ 3 

.... —^S « given that-y = 1 when or = 0. 

step—length /z ~ 0.T; we have- — 

= ?0 = 1 .' ~^f{*ry)^X +5F 3 . 


|fc&C3r,==p^.5^ 


^=(0.1X0+1) - 0.1, 


&.=&&* 0 +j,r 0 + k ^ 


=mmm+ U025) 
=@-1)11,1525) 
=$11525,./ 
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; h 


n 


Numerical Solutions of Ordinary Differential Equations 


h k 

/{ x Q + ^ ,j/ 0 + ^ 


= (0.1)(0.05 +1.1185) 

= (0.1X1.1685) 

= 0.11685 

k 4 = ^/(^o +h 'y 0 +k 3 ) 

= (0.1X0,01 +1.2474) 

= (0.1 )(1.3474) '/ ' 

= 0.13474, 

4y =!-(*, +2fc 3 + fc 4 ) 


= 2 (O.l I- 2(0.11525) + 2(0.11685) + 0.13 4 7 4) 

= - (0.6991)= 0.1165. 

6 V 


We get 


T-- 


y 1 = y 0 +z\y = 1 +0.1165 
/.jr(0.1)= 1.1165. 

For (he second step, we"have 
x 0 = 0.1,y Q = 1.1165, 
k, = (0.1X0.1 + 1.2466) = 0.1347, 

. k 2 = (0.1X0.15 +*T.4014) = (0.1)0-55S14) 
= 0.1551, 

k 3 = (ai)(0.15 + 1.4259) — (0.1 )(1.5759) 
z =0.1576, 

k 4 = (0.1)(0.2 + 1.6233) = (0.1X1-8233) 

= 0.1823, 

Ay =-(o:9.424) =0.1571, 

6 

y(0.2) = 1.1165+0.1571 = 1.2736' 
y(0.1) = 1.1165 andy(0.2) = 1.2736. 
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w ^ 1. .Solve the'equation ” ^ * ~~3 * $(&) = 1 for x =#.2 and-* fc>3 decimal-pfef^es 

by Runge-^iKutta^iouith o?derme£ha$.' 

2. Use theRunge—lCuUa method to 8tl a^d*=&2for the. eqrasa&ro** 

P>< - r ' - 

dy 

— = *+y -y(0>=l. 

dx - , 

- ^ 

. tfy . .. *V . / - . . • x 

& For the equation ~ ^ ~ , y£9) — 1. ifed y.aS the points. wn&h the 

given step — lengths. ^ 

4 Use Runge — Kutta ntefiSd; toisolVe y 5 ~ xyf&J-x = I.4y x = l,y ~ 2 {by taking; 

step — length Jr=0.2). 


dy y —2x 


use Range — Kutia method thffady at-*— &I>0fc2,03and (44, given 


d* y~-t-x 
ihaty = T when:r = 0- * 

5. Use Runge — Ku tta method to obtain y wftesx?=J-4g5V'ers-&^a9y = 1.2 when x = l and 
d$ ^ 2 

y satisfies the equation "— - ^ -r y . 

dx 


3k &-851,0.780 2.1.1103; 1/2928 
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a * 
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* G.l 

ox 1 

* i .0665242 
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oiv 

1.16?2m 


0.4 

. 0A 

1.4762 


2.994S5066 5: y(0.1) = 1.0&7* 
a *Q.I) = 1.7271 


y&42r> == y@0i= £23& # y®.4) = 1.2544 
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